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Abstract: Observations with the Fermi Large Area Telescope (LAT) indicate an ex¬ 
cess in gamma rays originating from the center of our Galaxy. A possible explanation for 
this excess is the annihilation of Dark Matter particles. We have investigated the anni¬ 
hilation of neutralinos as Dark Matter candidates within the phenomenological Minimal 
Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used 
to search for solutions within the pMSSM. We found solutions that are consistent with 
astroparticle physics and collider experiments, and provide a fit to the energy spectrum 
of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with 
a mass in the range 84 — 92 GeV or 87 — 97 GeV annihilating into W bosons. A third 
solutions is found for a neutralino of mass 174 — 187 GeV annihilating into top quarks. 
The best solutions yield a Dark Matter relic density 0.06 < Ulh? < 0.13. These pMSSM 
solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If 
the pMSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be 
discovered soon. 
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1 Introduction 

Observations of our Galaxy and other individual galaxies [1, 2], clusters of galaxies, gravi¬ 
tational leasing by clusters [3] as well as the detailed properties of the Cosmic Microwave 
Background [4] all infer that the mass density in the Universe (excluding the vacuum den¬ 
sity) is dominated by an unseen component; Dark Matter (DM). Current observational 
evidence, as well as considerations of standard Big Bang primordial nucleosynthesis, rule 
out that this unseen component is baryonic in nature, such as a large population of black 
holes or brown dwarfs [5]. 
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The most likely explanation therefore is that DM consists of a neutral, very weakly inter¬ 
acting particle outside the Standard Model of particle physics, with the currently leading 
hypothesis being Weakly Interacting Massive Particles (WIMPs) [6-9]. If this particle is 
a thermal relic, with a mass on the weak scale ~ 100 GeV, the velocity-weighted 
cross section should be of the order (av) ~ (2 — 5) x 10“^® cm^ s“^ [10, 11] in or¬ 
der to produce a DM density corresponding to — 0.12 as required by observa¬ 

tions (e.g. [4]). Here Ddm is the dark matter density in units of the critical density and 
h = Ho/{100 km/s per Mpc) ~ 0.68 with Hq the Hubble constant. 

Large-scale simulations of galaxy formation in the context of a flat ACDM cosmology all 
predict extensive, centrally concentrated, dark matter halos around galaxies such as our 
own [12, 13]. This implies that the strongest possible indirect DM signal should come from 
the Galactic Center (GC), in particular in the form of gamma rays from DM annihilation 
(for a recent review see [14]). Gamma rays with photon energies below 100 GeV are 
not attenuated or deflected during their flight over ~ 8.5 kpc from the GC, unlike other 
observable decay products [15]. 

Observations of the GC region with the Fermi-LAT satellite show a gamma ray excess 
for photon energies that peak in the range 1 GeV < < 5 GeV after a careful (and 

non-trivial) subtraction of the diffuse emission from known astrophysical sources [16-27]. 
These include gamma rays due to bremsstrahlung and from the decay of neutral pions 
produced by cosmic rays in the interstellar gas around the GC. The GC excess extends 
well away (> 10°) from the Galactic plane, as expected for a DM signal [24, 28, 29]. 
Therefore, even though a scenario where the GC excess is caused by conventional sources 
(e.g. unresolved point sources [30-34] or burst events associated with the 2 x 10® Mq central 
black hole [35, 36]) can not be completely excluded, a DM origin seems not unlikely. 
Other indirect searches with positrons [37, 38], anti-protons [39-45] or dwarf spheroidal 
observations [46-49] become increasingly sensitive to the required cross sections. 

There have been already a large number of attempts to explain the excess in a plethora 
of particle physics theories/models [50-100], including supersymmetric (SUSY) [101-114] 
scenarios [115]. Particular emphasis has been put in SUSY realizations beyond the minimal 
supersymmetric standard model (MSSM) [82, 91, 116, 117]. The reason is that in the 
MSSM, the required neutralino annihilation rate to the two golden channels, namely to 
T~^T~ and to bb with neutralino masses of ~ 10 GeV and ~ 30 GeV respectively (as found 
in most earlier analyses of the excess spectrum) is in tension with LEP or LHC bounds on 
sfermion masses. 

However, recently it has been shown that accounting for systematic uncertainties in the 
modeling of astrophysical backgrounds [118] opens up the possibility that the annihilation 
to other final states can fit the excess relatively well, even for DM masses as high as ~ 126 
GeV (in the case of h®h® final states) [115, 119]. This renews the interest in the question 
of whether the GeV excess can already be accommodated in the MSSM. 

In this paper we show how the MSSM offers explanations of the GC excess and how these 
scenarios ar going to be proved in the run H of the LHC and in the near future with the 
ton-scale DM direct detection experiments and in a complementary way by IceCube with 
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the 86-strings configuration. 

The paper is organized as follows. We describe the uncertainties involved in the GC excess 
in Section 2. In Section 3 we introduce our theoretical model and the methodology used for 
its exploration. Section 4 is devoted to present our results and Section 5 for our conclusions. 
Uncertainties in modelling the photon excess spectrum are discussed in the appendix. 

2 Galactic center observations in light of foreground systematics 

The observed gamma-ray flux from DM annihilation per unit solid angle at some photon 
energy is given by 



( 2 . 1 ) 


where the integral is along the line of sight (LOS) at an angle 9 towards GC, {av) is the 


(relative) velocity weighted averaged annihilation cross section, tudm denotes the DM mass, 
and dN/dE is the photon spectrum per annihilation. The flux is sensitive to uncertainties 
in the distribution in the radial DM density profile, /9 dm(?’), as function of galactocentric 
distance r. Dark matter-only simulations of large-scale galaxy formation can in principle 
resolve the central ~ 1-2 kpc of DM halo (e.g. [120]). However, for our Galaxy, DM 
dominates the dynamical estimates for the total (baryonic -|- DM) enclosed mass, M(< 
r) ~ only beyond a galactocentric distance of 20 kpc, as can be obtained from 

galaxy rotation curves Uot(^)- This renders the inner DM density profile rather uncertain, 
see for instance [121]. 

It is quite common to adopt a generalized Navarro, Frenk & White (NFW) profile [12], 
with /?dm(^’) oc (r -|- rs)“~^, with a = 1 for the original NFW profile. The radius is 
usually taken to be around 20 kpc, which implies ^ close to the GC. 

The main uncertainties are twofold: (1) Infall of baryonic gas towards the GC in the late 
stages of galaxy formation initially steepens the DM density profile, increasing a, while 
mass loss due to supernova-driven winds from the first generation(s) of massive stars in the 
Galactic Bulge can flatten it. The net effect is difficult to determine in general, but recent 
simulations that combine DM with hydrodynamics for the baryonic content [122] show a 
flattening of the density profile for Mily Way like spiral galaxies (2) The normalization of 
the DM density distribution is difficult to determine. It is usually parametrized by the 
DM density at the galactocentric distance of the Sun, pdm(?'©)• Global determinations 
and local determinations in the Solar neighborhood yield values in the range /Odm(?'©) — 
0.2 —0.5GeV/cm^. The main uncertainties in global determinations stem from modeling of 
the shape of the halo, while local determinations suffer from uncertainties in the baryonic 
surface density of the Galactic disk and/or the local stellar kinematics [123, 124]. 

The consequence for predictions of the flux of the GC excess is that, with particle physics 
parameters fixed, the uncertainty in the predicted absolute flux level exceeds a factor of 
a few for realistic parameters. Throughout, we will adopt the estimates of the J-value 
uncertainty as discussed in [119]. There, the uncertainty of the signal flux at 5 degree dis¬ 
tance from the Galactic center was estimated by scanning over a large range of generalized 
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NFW profiles that are consistent at the 95% CL with the microlensing and rotation curve 
constraints from [125]. The corresponding J-value uncertainty is (very conservatively, since 
additional constraints from the slope of the profile in the inner 5 degree are not taken into 
account) a factor of ~ 5 in both directions. 

The existence of a spectrally broad and spatially extended “excess” emission (“Fermi GeV 
excess”) above conventional convection-reacceleration models for the diffuse gamma-ray 
emission is by now well established. One of the possible explanations that can explain the 
properties of this emission surprisingly well is the emission from the annihilation of DM 
particles. 

In order to search for corroborating evidence for the dark matter interpretation of the 
excess, it is important to estimate the uncertainties of its spectral properties conservatively. 
We adopt here the results from [118], where the excess emission was studied at latitudes 
above 2 degree. This region is very sensitive to a dark matter signal, but avoids the much 
more complicated Galactic center region. The corresponding likelihood function will be 
discussed below in Section 3. 

The MSSM is still the most promising framework for WIMP dark matter models. However, 
as we will show, it is not completely trivial to find valid model points which provide a 
spot-on description of the spectrum of the GeV excess. However, in order to not dismiss 
possible collider signatures that would serve as corroborating evidence for a dark matter 
interpretation, we will allow below for additional uncorrelated systematics that might affect 
the spectrum and discuss additional uncertainties e.g. coming from the predictions of the 
photon energy spectrum from dark matter annihilation, as discussed below. In the case 
that the DM origin of the GeV excess is supported by other experiments, these additional 
uncertainties require further study. 

3 Analysis setup 
3.1 The Model 

The MSSM has 105 Langrangian parameters, including complex phases. One can reduce 
this number to 22 by using phenomenological constraints, which defines the so-called phe¬ 
nomenological MSSM (pMSSM) [126]. In this scheme, one assumes that: (i) All the soft 
SUSY-breaking parameters are real, therefore the only source of CP-violation is the CKM 
matrix, (ii) The matrices of the sfermion masses and the trilinear couplings are diago¬ 
nal, in order to avoid FCNCs at the tree-level, (hi) First and second sfermion generation 
universality to avoid severe constraints, for instance, from mixing. This number 

can be further simplified to 19 parameters (we will refer to this here as pMSSM) and still 
capture the phenomenology of the 22-parameter model. 

The 19 remaining parameters are 10 sfermion masses,^ 3 gaugino masses Mi^ 2,3 j the ratio 

^The corresponding sfermion labels are Q\, Qs, L\, L 3 , wi, di, U 3 , d^, ei and 63 . Here 1 indicates the 
light-flavoured mass-degenerate 1st and 2nd generation sfermions and 3 the heavy-flavoured 3rd generation. 
The labels Q and L refer to the superpartners of the left-handed fermionic SU{2) doublets, whereas the 
other labels refer to the superpartners of the right-handed fermionic SU{2) singlets. 


~ 4 - 



of the Higgs vacuum expectation values tan/3, the Higgsino mixing parameter /i, the mass 
ruA of the CP-odd Higgs-boson and 3 trilinear scalar couplings Ab^t,T- 
In this scenario, in principle, there are five arbitrary phases embedded in the parameters 
Mi{i = 1, 2, 3), fj, and the one corresponding to the trilinear couplings provided we assume 
that the trilinear matrices are flavour diagonal. However one may perform a U{1)r rotation 
on the gaugino fields to remove one of the phases of Mj. We choose the phase of M 3 to 
be zero. Note that this U{1)r transformation affects neither the phase of the trilinear 
couplings, since the Yukawa matrices being real fixes the phases of the same fields that 
couple to the trilinear couplings, nor the phase of /i. Therefore in the CP-conservation case 
Ml, M 2 , /r and the trilinear couplings can be chosen both positive and negative. 

3.2 Generation and pre-selection of pMSSM model-sets 

For our exploration of the pMSSM we use SUSPECT [126] as spectrum generator. Dark- 
SUSY 5.1.1 [127, 128] is used for the computation of the photon fluxes and MicrOMEGAs 

3. 6 .9.2 [129, 130] to compute the abundance of dark matter and and cr®5p. 

For the hadronic matrix elements /t„, /t^ which enter into the evaluation of 

the spin-independent elastic scattering cross section we adopt the central values presented 
in Ref. [131]: = 0.0457, = 0.0457. For the strange content of the nucleon we use 

recently determined average of various lattice QCD (LQCD) calculations = 0.043 [132]. 
The spin-dependent neutralino-proton scattering cross-section depends on the contribution 
of the light quarks to the total proton spin A„, and A^. For these quantities, we 
use results from a LQCD computation presented in [133], namely A^ = 0.787 ± 0.158, 
Arf = —0.319 ± 0.066, As = —0.02 ± 0.011 [133] and leave them vary in the la range. We 
will explain why we adopt this approach later. 

Following [134], we assume that the ratio of the local neutralino and total dark matter 
densities is equal to that for the cosmic abundances, thus we adopt the scaling Ansatz 

? =(3.1) 

For Hdm we adopt the central value measured by Planck, Hdm = 0.1186 [135]. The photon 
fluxes are rescaled with when the predicted value is below 0.0938 which encompasses the 
2a level uncertainties both in the theoretical prediction and the value inferred by Planck 
added in quadrature. This allows multi-component Dark Matter. 

We select only models with a neutralino as lightest SUSY particle (LSP). From SUSY 
searchers at colliders we impose the LEP limits on the mass of the lightest chargino. Namely 
103.5 GeV [136]. The Higgs mass has been precisely determined by ATLAS and GMS to 
be 125.4 (ATLAS [137]) and 125.0 GeV (CMS [138]) with uncertainties of 0.3 — 0.4 GeV 
for each experiment. On top we account for a theoretical error of 3 GeV [139] in its 
determination and select models with a lightest Higgs boson within the range: 

122 GeV < rriho < 128 GeV . (3.2) 

From the dark matter point of view we in addition demand the following constraints: 

• Upper limits from the LUX experiment on the spin-independent cross section [140]. 
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• Upper limits from the IceCube experiment with the 79 string configuration on the 
spin-dependent cross section [141], assuming that neutralinos annihilate exclusively 
to W~^W~ pairs. 

In the parameter scan it was required that solutions need to have Ma > 800 GeV or 
5 < tan(/3) < 0.075 • Ma — 16.17 to ensure that they are not excluded by searches for heavy 
Higgs bosons. 

3.3 Parameter scan 

In a first iteration the pMSSM parameter space was randomly sampled with > 10® param¬ 
eter points from a flat prior. All possible DM annihilation channels have been compared 
to the measured Fermi photon flux in two energy bins around 1 and 5 GeV. All mass 
parameters were sampled between —4 TeV and 4 TeV. 

In an iterative procedure the best ht points of the first iteration were used as seeds to sample 
new model parameter ranges centered around the seed points and with multi-dimension 
Gaussian distribution as widths. The ranges of some parameters were reduced: 100 GeV 
to 1 TeV and —1000 GeV to —100 GeV for Mi and M 2 , 100 GeV to 1000 GeV for n and 
tan/3 between 1 and 60. The iterative sampling procedure was repeated several times, until 
a reasonable annihilation process was found. The process was found to be XiXi W^W~ 
for our first and second solution and XiXi for the third solution. The main annihilation 
diagram is the t-channel exchange of a xf (or the t-channel exchange of a stop quark). 

In the final iterations 11 of the 19 parameters have been set high enough to be non-relevant 
(4 TeV). The final set of parameters influence electroweakinos, the Higgs mass and the spin- 
independent cross section. The final set of parameters was: 

Ml, M 2 , n, tan^S, Ma, ds, Q 3 , At- 

3.4 Galactic Center excess region 

For all model points DarkSUSY was used to derive the photon spectrum dN/dE of the 
annihilaton process, which was then compared to the spectrum of the GeV excess emission. 
We adopt the x^ dehnition from [118], which takes into account correlated uncertainties 
from the subtraction of Galactic diffuse gamma-ray backgrounds. However, in addition to 
the astrophysical uncertainties in the measured spectrum as discussed in [118], we allow 
for an additional 10% uncorrelated uncertainty in the predicted spectrum, as motivated in 
Appendix A. 

We use the following definition 

X^ = '^{di - mi){T.ij)-^{dj - rrij) , 
ij 

where i and j are the energy bin numbers running from 1 to 24, di and rrii is the Fermi 
and model flux, respectively, and Sjj is the covariance matrix that incorporates all relevant 
statistical and systematic uncertainties when modeling the GeV excess flux. As mentioned 
above, we will allow for an additional uncorrelated systematic uncertainty of the level 
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Figure 1: Photon excess spectrum as extracted in Ref. [118] from the Fermi data from 
the inner Galaxy, compared with the model calculations with the lowest Xio figure, 
p-value= 0.3 with x^q) and the model with the lowest Xo (I'ight figure, p-value= 0.025 with 
Xq), for WW solution 1. Note that besides the statistical errors, which are shown as error 
bars, there are two kinds of systematics which affect the observed photon spectrum (shown 
as gray dots): Firstly, there are uncertainties from the removal of astrophysical foregrounds 
(shown by the gray boxes; mostly inverse Compton and tt® emission, see Ref. [118] for 
details). These uncertainties are strongly correlated and can lead in general to an overall 
shift of all data points up or down, as illustrated by the black dots. Secondly, there are 
particle physics uncertainties in the predicted photon spectrum, which we conservatively 
assume to be at the 10% level (green band in left panel, only affecting Xio)- Details are 
discussed in Appendix A. 


of (Tg = 10%, which is incorporated in the covariance matrix from [118] by substituting 
Tiij —)■ Sjj + dijdfa^. Photon generation via hadronic or top decays is mainly caused 
by Quantum Chromo Dynamic processes which are described with semi-empirical models 
with many parameters. Also the uncertainties in the photon energy scale can change the 
shape in the modelling of the photon excess spectrum (see Appendix A). 

In the following Xo denotes as = 0% and Xio denotes ag = 10%. Some distributions are 
shown with both definitions to illustrate the effect of including uncorrelated systematic 
uncertainties in the predicted photon spectrum. 

4 Results 

4.1 The galactic center excess 

In our exploration of the pMSSM parameter space we find that requiring a Xio < 
(corresponding to a p-value > 0.02) implies the following three pMSSM parameter ranges: 

4.1.1 WW solution 1: Bino-Higgsino neutralino 

In this type of solution, the neutralinos annihilate mostly exclusively to W~^W~ pairs. 
Only a small fraction annihilate to W~^W~/bb. The reason is that even being away of the 
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Figure 2: The neutralino mass as a function of Mi, M 2 and /r. is shown as colour 
code. 


A-funnel region the neutralino coupling to pseudoscalars is enhanced due to their bino- 
higgsino nature and therefore their annihilation to pairs of b-quarks. 

This solution provides a good (and in our scan the best) fit to the Galactic center photon 
spectrum as measured by Fermi. This is partly due to the fact that we, in contrast to 
previous studies, allow for an additional 10% uncorrelated uncertainty on the predicted 
photon energy spectrum, as discussed and motivated in Appendix A. The best fit points 
have Xio ~ 27 (p-value ss 0.3) with the best-fit normalization of the Xo'^^ and a Xio ~ 24 
(p-value=0.45) with the best-ht normalization of the Xio'^f (here we take 10% uncertainties 
in the predicted spectrum into account in the ht, see above). The best Xo found to 
be ~ 39.5. Figure 1 compares the photon spectrum as measured by Fermi with the model 
calculations with the lowest xfo Xo- 

The properties of these models are shown in Figure 2 (tagged as WW(1). The best solutions 
correspond to: 

-103 GeV < Ml < -119 GeV, 

240 < M 2 < 660 GeV, 

108 GeV < /i < 142 GeV, 

8 < tan < 50. 
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Figure 3: Q.h^ as a function of the mass of the DM candidate, is shown as colour code. 
Both definitions are shown. 


It can be notice that the bino mass Mi and the higgsino mass /U are very strictly constrained 
leading to a precise forecast for DM direct/indirect detection and LHC physics. 

The composition of the lightest neutralino is ~ 50% bino and ~ 50% higgsino and the mass 
is in the range ~ 84 — 92 GeV. 

Figure 3 shows that all points tagged as WW(1) with Xio < 35 correspond to Q-h? in the 
range ~ 0.07 — 0.125. Recall that this constraint was not used in the fit procedure. We 
consider the outcome as remarkable since can vary between ~ 10 “^ and ~ 10 ^ within 
pMSSM models. 

In terms of contraints coming from electroweakino searches at the LHC M 2 is less tightly 
constraint and ranges between about 300 — 900 GeV. If M 2 is smaller than about 170 — 250 
GeV, the corresponding neutralino (the X 4 ) decays to Z and Xi- This little part of the valid 
parameter region is excluded by LHC chargino-neutralino searches already. If M 2 > 250 
GeV the X 4 decays into charginos, Z and Higgs bosons. This region is not much constrained 
at the LHC so far. LHC signatures are further discussed in the next section. 

Finally, Fig. 5 shows that points consistent with this solution have a pseudoscalar mass 
^ 350 GeV, therefore the points that fit well the GC excess lie to the SUSY decou¬ 
pling regime in which the lightest Higgs is Standard Model like, thus consistent with LHC 
measurements of the Higgs properties. 

4.1.2 WW solution 2: Bino-Wino-Higgsino neutralino 

As in the case above, in this type of solution the neutralinos annihilate mostly exclusively 
to W+W- pairs. The following parameter range yields p-values between 0.02 and 0.15: 

91 GeV < Ml < 101 GeV, 

102 GeV < Ms < 127 GeV, 

156 GeV < /i < 507 GeV, 

5 < tan/3 < 12 
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The composition of the neutralino is dominant bino (~ 90%) with a ~ 6 % of wino and 
a ~ 4% of higgsino whereas the mass is in the range ~ 86.6 — 97 GeV. Figure 3 shows 
Vth? as a function of the mass of the DM candidate (points tagged as WW(2) ) with the 
corresponding The best ht points have 0.05 < < 0.15 consistent with Planck. 

The LHC sensitivity to this scenario is similar to the Bino-Higgsino case since the only 
difference is that in this case the neutralinos Xg 4 are heavier than the others. Figure 5 
shows, as in the Bino-Higgsino solutions, that the lightest Higgs is “Standard Model like”. 

4.1.3 Top pair solution 

The third solution yields mostly neutralino annihilation into a pair of top quarks via the 
t-channel exchange of a right-handed stop quark. The neutralino is mostly Bino ~ 99% 
and in this case the chirality suppression in the annihilation cross section that affects to 
the other fermion final states does not apply here. 

As displayed in Figure 2 the solutions (tagged as tt) have a maximum p-value of 0.1. The 
best solutions imply the following pMSSM parameter range: 

171 GeV < |Mi| < 189GeV, 

190 GeV < IM 2 I < 1550GeV, 
p > 250 GeV, 
tan /3 > 5 

The neutralino mass is about the kinematical threshold ~ 174 — 187 GeV and the 
right-handed stops have a mass of ~ 200 — 250 GeV whereas the left-handed are heavy 
with a mass ~ 2600 — 3700 GeV to fulfill the Higgs mass constraint. 

In this case, as it can be seen Figure 3, all points tagged as tt cover a wider range than in 
the previous solutions for (~ 0.066 — 0 . 22 ). 

The right-handed stops decay to the lighter chargino and a bottom quark. The chargino is 
close in mass with the lightest neutralino (A ~ 50 GeV) leading to a hardly visible signal. 
Therefore this scenario evades current LHC constraints from stop searches. 

As above. Figure 5 shows that the pseudoscalar mass ruA > 500 GeV, therefore the lightest 
Higgs is Standard Model like. Figure 8 summarizes the third generation parameters found 
in the different solutions. The scan localizes very small volume elements of the parameter 
space. 

4.2 Implications for DM direct and indirect experiments 

Dwarf spheroidal galaxies New recent observations of dwarf spheroidal galaxies with 
the Fermi Large Area Telescope provide by now the most stringent and robust con¬ 
straints on the velocity-averaged annihilation cross-section [142]. These limits are usu¬ 
ally considered have to be taken into account when interpreting the emission seen from 
the Galactic center in terms of dark matter annihilation. The for us most relevant final 
states are W~^W~; for a dark matter mass around 80-90 GeV, current upper limits are 
{av) < 2.6 X lO-^^cm^s-i [142]. 
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Figure 4: The velocity averaged annihilation cross section, {av) as a function of the mass 
of the DM candidate. Xio is shown as colour code. We also show the 95%CL upper limits 
obtained from a combined observation of dwarf spheroidal galaxies in Ref. [142]. 


As can be seen from Fig. 4, this constraint is fulfilled by the models considered in this 
work. In fact, all interpretations presented in this paper require a relatively large J-value 
at the Galactic center, which implies annihilation cross-sections that are smaller than the 
thermal value. Hence, although dwarf spheroidal observations could potentially confirm a 
dark matter interpration of the GC excess in the future, they cannot currently be used to 
rule out an interpretation in terms of the MSSM. 

Spin-dependent and spin-independent cross sections Within the MSSM the dom¬ 
inant contribution to the spin-independent (SI) cross-section amplitude, when squarks are 
heavy, is the exchange of the two neutral Higgs bosons. The SI cross-section for H/h ex¬ 
change is oc |(A^i 2 — A"!! tan 0 u,)p|A';^ 3 /]^ 4 |, where 9w is the electroweak mixing angle, Nu 
represent the neutralino composition. 

With regard to the spin-dependent (SD) cross-section, the dominant contribution corre¬ 
sponds to the exchange of a Z boson. Since the bino and wino are both SU(2) singlets, 
they do not couple to the Z boson, and therefore SD cross-section is largely determined by 
the higgsino content of the neutralino. The Z exchange contribution (and hence the SD 
cross-section) is proportional to the higgsino asymmetry (| Aisp — |Ai 4 p)^. The asymmetry 
is maximized when either the binos and higgsinos or winos and higgsinos are close in mass. 

4.2.1 WW solution 1: Bino-Higgsino neutralino 

In solutions of the bino-higgsino type one expects large SI cross-sections as explained above. 
In fact, the lightest Higgs contribution is effectively fixed and pushes the SI cross-section to 
values that are in conflict with LUX bounds, therefore cancellations with the heavy Higgs 
are required. It is well known that these cancellations arise in non-universal models [143]. 
The degree of cancellation spans the SI cross-section down to ~ 10“^^ pb. Those cross 
sections are going to be probed by ton-scale experiments as Xenon. 
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This can be seen in the left panel of Figure 6 (points tagged as WW(1)) where we show the 
(c^Lp, rriy^) plane with the current 90% exclusion limits from the LUX collaboration. The 
result is rescaled with the scaling Ansatz of Eq. (3.1) to account for the fact that the local 
matter density might be far less than the usually assumed value local p = 0.3 GeV/cm^. 
In the right panel of 6 we display the (cr^^pj rux) plane with the current 90% exclusion 
limits from the IceCube collaboration with the 79 strings configuration assuming that 
the neutralinos annihilate exclusively to W~^W~ [141]. Here the SD cross section is not 
rescaled since the IceCube detection depends on whether the Sun has equilibrated its core 
abundance between capture rate and annihilation rate. Typically for the Sun, equilibration 
is reached in our points. 

Since the higgsino asymmetry is sizable in this scenario, the SD cross-sections are large 
and close to the current limits imposed by IceCube. Actually, the model becomes tightly 
constrained and one has to allow, at least, up la deviation of the central values for the 
hadronic nucleon matrix elements for SD WIMP nucleon cross sections estimated using 
LQCD. It is interesting to notice that all the currently found points are within the reach 
of IceCube with the 86 strings configuration. Therefore this phase space is going to be 
probed in a near future. 

4.2.2 WW solution 2: Bino-Wino-Higgsino neutralino 

These type of solutions are expected to follow a similar pattern to the Bino-Higgsino 
scenario. Specially in terms of the SI cross section. This is verihed in the left-panel of 
Figure 6 (points tagged as WW(2)) from where one can infer that ton-scale experiments 
will probe a sizable fraction of the parameter space consistent with this scenario. 

The fact that the Higgsino composition is reduced alleviates the tension in the SD cross 
section with respect to the current bounds set by IceCube as it can be seen in the right-panel 
of Figure 6 (points tagged as WW(2)). Indeed we find that all our points are well below the 
current IceCube limits even taking central values for the hadronic nucleon matrix elements 
for the SD WIMP nucleon cross sections estimated using LQCD. In terms of prospects 
most of the points are out of the IceCube reach. 
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Figure 6: a si (left-panel) and asD (right-panel) as a function of the mass of the DM 
candidate, is shown as colour code. 


4.2.3 Top pair solution 

With regard to DM detection, points lying to this scenario are expected to have different 
features with respect to the previous type of solution because the neutralino is mostly bino 
~ 99%. It leads to a lower prediction for both the SI and SD cross sections as it can be 
seen in both panels of Figure 6 (points tagged as tt). The most evident differences arise 
in the SD cross section which now expands down to values of r-u 10 pb. Clearly this 
scenario is not going to be fully proved for experiments sensitive, both, to SI and SD cross 
sections. Despite this, experiments sensitive to the SI cross sections as Xenon 1-ton will 
probe some fraction of the parameter space consistent with this scenario. 

4.3 Implications for LHC searches 

4.3.1 WW solution 1: Bino-Higgsino neutralino 

Since the neutralino and chargino mixing matrix parameters are highly constrained in the 
allowed parameter region the production rates and decays of all neutralinos and charginos 
are constrained. 

Neutralino X 123 Higgsinos and Binos, the x^ is a Higgsino. All these electroweakinos 
have very similar masses. The decay of the ^2 3 x^ to the LSP will not lead to high 
energetic signals. Consequently the production of the 3 light Neutralinos and the light 
Chargino will not be visible at LHC in neutralino-chargino searches. 

We see a few interesting LHC signals: 

Chargino-t-Neutralino production. The only signal visible in electroweakino searches 
at the LHC could be x^x^ production with the subsequent decays of X 4 to Zxi, Higgs-t-x)* 
and W-|-xf ■ Higgs production in this scenario is discussed in [144]. 

Monojets. Since the lightest 3 neutralinos have a similar mass and a Higgsino component 
they can be pair produced via s-channel Z production. In addition the xf can be produced. 
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The combined cross sections is enhanced compared to XiXi alone. This might lead to a 
signal in monojet events for the upcoming LHC data. 

Searches for squarks and gluinos. Finally searches for squarks and gluinos can be 
conducted in our scenario. If Mi, M2,/i, tan/3 are fixed, the decays of squarks and gluinos 
is well determined yielding specific signatures. Especially right-handed squarks will likely 
decay via the heavy Winos leading again to Z and Higgs signals. 

4.3.2 WW solution 2: Bino-Wino-Higgsino neutralino 

For this solution Neutralino Xi 2 mostly Bino-Wino and have a masses of ~ 88 and 
~ 106 GeV. xf is mostly Wino and has a mass of ~ 105 GeV. On the other side there are 
3 heavier states {X 34 and xf) with a mass of around 400 GeV. 

The LHC signaturs are similar to solution 1; 

Chargino-|-Neutralino production. The three heavier states will be visible in the 
searches for chargino-neutralino production. Again the heavy neutralinos will decay into 
ZXi, Higgs-Fy? and W+xf- 

Monojets. Since the lightest 2 neutralinos and the lightest chargino have a similar mass 
and a Higgsino component they will be visible in monojet production. The cross section 
will be small compared to solution 1 and the signal will be harder to detect. 

Searches for squarks and gluinos. For squark and gluino searches the conclusion is 
similar to solution 1. 

4.3.3 Top pair solution 

Interesting is that also our third solution seems also not excluded by run-1 LHC searches. 
The neutralino x^ 2 again mostly Bino-Wino and have a masses of ~ 170 and ~ 225 
GeV. xf is mostly Wino and has a mass of ~ 225 GeV. Again we have 3 heavy (dominantly 
higgsino) states {X 3 4 and xf) with a mass of around 850 GeV. 

The solution will lead to the following signatures for run 2: 

Chargino-|-Neutralino production. The light neutralino states are again quite com¬ 
pressed and might only be visible with a very soft lepton search. 

Monojets. The compressed light neutralinos and chargino have masses of ~ 170 GeV 
which reduces the cross sections for monojet searches compared to the WW scenarios 
discussed above. 

Search for stops pair production. The stop mass is ~ 230 GeV. The stop decays 
100% to xf and a b-jet. The xf has a mass difference of 50 GeV with the Xv This 
signal should be visible with dedicated stop searches in the upcoming run-2 data. 
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the mass of the DM candidate, is shown as colour code. 


4.4 Implications for flavour observables 

Finally in this section we discuss the implications for flavour physics. In Figure 7 we show 
on the left-panel the BR{Bs —)• and on the right one the BR[B —)• X^y) versus the 

neutralino mass. 

Accounting for both parametric and theoretical uncertainties in both observables and 
adding them in quadrature to the experimental ones implies that the allowed range at 
2(7 level is [145]: 

1.39 X 10"^ < BR(Bs < 4.49 x 10"^ 

2.76 X lO""^ < BR{B X^y) < 4.34 x 10"^. 

Let us first discuss the BR{Bs —)■ In the left-panel of Figure 7 one can see that all 

points corresponding to, both, the Bino-Higgsino and Bino-Wino-Higgsino neutralino type 
of solutions are within the range above. This is quite remarkable since we have not used 
this observable as constrained in our scan. 

In the top pair type of solution the conclusion is broadly the same with the exception 
of a few points which are ruled out. Those correspond to tan /3 > 40 where new physics 
contributions are sizable in the minimal flavour violation scenario [146]. In particular, when 
stop quarks are relatively light. This is precisely which makes the distinction between the 
Bino-Higgsino, Bino-Wino-Higgsino neutralino and top pair type of solutions as it has been 
already pointed out. 

In the BR[B —)• Xs'y) case, the results are shown in the right-panel of Figure 7. Here a 
fraction of the points belonging to the Bino-Higgsino neutralino solution are ruled out by 
current experimental bounds whereas most of points corresponding to both the Bino-Wino- 
Higgsino and top pair solutions are allowed. The largest values correspond to relatively 
large tan /3 values together with the fact that the lightest chargino is Higgsino like and the 
interference with the Standard Model contribution is positive since sgn{fj,At) > 0 [147]. 
Again it is worth stressing that most of the solutions are allowed without imposing this 
constraint in our scan. 
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5 Discussion 


We have systematically searched for Dark Matter annihilation processes to explain the 
excess found in the photon spectrum of the Fermi-LAT satellite. We found three solutions 
where the excess is explained by the annihilation of neutralinos with a mass around 84 — 
92 GeV, 86 - 97 GeV or 174 - 187 GeV. 

These solutions yield the following interesting features: 

• The neutralino of our first and second solutution is a Bino-Higgsino or a Bino-Wino- 
Higgsino mixture annihilating into W~^W~. We obtain a good fit to the Galactic 
center gamma-ray data by allowing for an additional (and reasonable) uncertainty of 
the predicted photon spectrum of 10%. The corresponding neutralino and chargino 
mixing parameters are well constrained for both solutions. 

• A third solution is found where a (dominantly Bino) neutralino annihilates into tt, 
which provides however smaller fit probability for the Galactic center data. 

• Since light electroweakinos are compressed, this sector is hard to test at the LHG, but 
might lead to a signal in monojet (or soft-lepton monojet) events in the upcoming 
LHC run. In addition the production of the heavy Wino (or mixed) states will be 
visible for most models. 

• Part of the spin-independent cross section can be probed by the upcoming ton-scale 
direct detection experiments. 

• All models points with a Bino-Higgsino neutralino have spin-dependent cross section 
which are well in the reach of the upcoming spin-dependent constraints provided e.g. 
by IceCube. 

• The best solutions yield values with 0.06 < fih^ < 0.13. This is a remarkable feature 
since fih^ varies for pMSSM solutions unconstrained by the Galactic center excess 
by about 10 orders of magnitude. 

If the MSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might 
be discovered soon. The solutions also exist in extensions of the MSSM with a similar stop 
and electroweakino sector. 
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A Uncertainties in the predicted photon spectrum 

We discuss here briefly sources for uncertainties in the predicted photon spectrum (see 
[148] for an earlier assessment), and leave a more detailed study to a future publication. 

Generation of the photon spectrum with Pythia. Dark Matter particles are not 
charged and cannot directly couple to photons. The Fermi-LAT excess spectrum can be 
described by Dark Matter (neutralino) annihilation to various SM particles (e.g. W^W~ 
in our models), which then decay further. The decay products can be quarks, which are 
influenced by the strong force. These quarks can further radiate gluons, which can split into 
further quarks. This is modelled within Monte Carlo event generators with semi-empiric 
models (e.g. so called Parton Showering). The quarks are then re-connected to colourless 
hadrons (again by models based on measurements of fragmentation functions). These 
hadrons decay and some have significant decay fractions to photons. The photon spectrum 
is given to a large amount by the momenta and multiplicity distributions of hadrons. By far 
most important are the decays of neutral particles (mainly vr*^), but photons are radiated 
at each moment in the chain. The spectrum of photons produced e.g. by decays has 
never been directly measured down to the energies relevant for the Fermi-LAT spectrum. 
The generation of a photon spectrum with Monte Carlo event generators has uncertainties 
stemming from the used model and the model parameters. Here we compare for the same 
generator and version (Pythia 8.1 [149]) various different hts of the model parameters (see 
also [150]). The photon spectra are shown in Figure 9 for the annihilation of neutrinos 
with an energy 85 GeV into W~^W~. Besides small effects stemming from the mass of the 
t-channel propagator the spectrum is identical with the annihilation of a DM particle with 
a mass 85 GeV into W~^W~. The differences range between 5-10% at low photon energies 
between 0.5-20 GeV and > 20% at larger energies. This uncertainty should be regarded 
as a lower limit, since no estimate was done to determine the parameter uncertainties via 
a full extrapolation of data uncertainties. Also no other models (as implemented e.g. in 
Herwig) have been considered. 
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Figure 9: Effect of a variation of the Pythia 8 tunes on the generated photon spectrum 
from vV —)• W^W~ with neutrino energies of 85 GeV. 


As discussed in the main text, the influence of such additional uncertainties is large: The 
best-MSSM fit has a p-value of 0.35 including a high-energy physics uncertainty of 10% 
and p-value of 0.03 without high-energy physics uncertainties. 

Variation of the photon energy scale. Another significant source of uncertainties is 
the uncertainty in the photon energy measurement of the Fermi LAT. The photon energy 
measurement has an uncertainty of 3 — 5% [151] measured in a range ~ 6 — 13 GeV. We 
assume a ± 1-sigma energy measurement uncertainty of ±5% for the unmeasured region 
3 — 5 GeV as reasonable. We determined the effect on the spectrum by changing the energy 
of each measured photon by -|-5% or —5% for all photon energies (and for comparison by 
± 10 %). 

Figure 10 shows the Pythia generated excess spectrum for neutrino annihilation into 
W+W- with a neutrino energy of 85 GeV. Nominally, the photon spectrum varies by 
± > 5% at energies of > 5 GeV. We conclude that such uncertainties need to be considered 
in the interpretations of the Fermi excess spectrum. However, we note that a photon energy 
rescaling does mostly affect the normalization, and not so much the shape of the spectrum. 
Since the change in the normalization is still much smaller than the uncertainties of the 
astrophysical J-value, the impact on the fit-quality is in fact not large: Only changing the 
fit-template from the nominal (no energy variation) to 5% up and 5% down changes Xq 
from 37.8 to 40.4 (up) or 35.3 (down). The p-value changes from 0.035 (nominal) to 0.02 
(up) or 0.065 (down). 
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Figure 10: Effect of a variation of the photon energy scale by ±5 — 10% on the generated 
photon spectrum from vv —?■ W^W~ with neutrino energies of 85 GeV. Note that the 
main effect is an overall change in the normalization (which has to be compared with the 
large uncertainties of the J-value) and a shift of the peak energy in log-space (which can 
mildly affect the quality of the fit). 
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